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Abstract 

We present numerical evidence that the hopping of electrons between chains 
described by the t — J model is coherent in the integrable cases ( J = and 
J = 2) and essentially incoherent otherwise. This effect is not related to 
the value of the exponent a, (which is restricted to the interval [0,1/8] when 
< J < 2), and we propose that enhanced coherence is characteristic of 
integrable systems. 
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After the proposal by Anderson that two-dimensional models of correlated electrons 
might behave as Luttinger liquids i.e. as their one-dimensional analogs, a lot of work 
has been devoted to studying the effect of a small transverse hopping t± between two chains. 
The first conclusion that has been obtained by several authors is that a small transverse 
hopping is a relevant perturbation as long as the exponent a describing the power-law sin- 
gularity of the momentum distribution in the isolated chains is smaller than 1, suggesting 
that a 2D system of coupled Hubbard chains flows toward a Fermi liquid fixed point, a 
being at most 1/8 in that case. More recently P], another point of view has been empha- 
sized according to which what really matters is not whether t± is relevant, but whether 
it induces a coherent hopping between the chains. A study of the short time dynamics of 
two coupled Luttinger liquids by Clarke, Strong and Anderson |^ suggests that coherence 
can be destroyed even if t±_ is relevant. Within their approximation, the hopping between 
chains becomes incoherent for a = 1/2, whereas t± remains relevant until a = 1. While this 
result opens new perspectives in the problem of coupled chains, this analysis suggests that 
the key parameter is still a, which somehow reflects the strength of the correlations within 
the chains. In a very recent paper about a related problem in other systems, Chakravarty 
and Rudnick ^ suggest that incoherence is actually a generic feature, and that coherence 
is limited to a very narrow range of parameters. The origin of this coherence is left as an 
open issue. 

In this Letter, we identify integrability as a very important factor of coherence between 
chains. More specifically, we show that the short time dynamics of the hopping of electrons 
between two chains is coherent if the model describing each chain is integrable, that this 
coherence is destroyed if one goes away from integrable points, and that this effect cannot 
be understood in terms of the exponent a. The system we have studied consists of two 
chains (so called 2 x L ladder) described by the t — J model and coupled by a perpendicular 
hopping t±. The Hamiltonian can be written 

a=l,2 \ i,<T i I i 
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where Pq is the Gutzwiller projection operator that excludes double occupancy on chain 
a. For small J/t ratios (say J/t < 0.5) an additional small transverse exchange coupling 
J± = J{^Y has also been included in order to improve the equivalence with the large- 
U Hubbard ladder in this parameter regime. In the isotropic regime {t±_ = t) transverse 
coherence can be established Q but hereafter we restrict ourselves to the small t± regime. 
Following Clarke et al, our analysis of coherence is based on the probability -P(t) that a 
system comes back after some time r to its initial state if one initially puts more particles 
on one chain than on the other. More precisely, -P(r) is defined by 

P{r) = \A{T)\' (2) 
A{t) = (^ole^^^-^^^l^o) (3) 

where {ipo) is the lowest energy [Eq) eigenstate of the system with t± set to zero and with 
AA^e more particles on one chain than on the other, while H is the full Hamiltonian. There 
are actually two such states because the excess particles could be on either chain. In the 
following, we have taken the symmetric combination of these states for numerical convenience 

i- 

In the case of non-interacting particles (referred to hereafter as U = 0) P{t) ex- 
hibits an oscillatory behavior characteristic of coherent transverse motion. More precisely, 
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P(r) = cos'^^'=(t_Lr) + i'^^" sin^^" (t^r) showing a characteristic period vr/ (4t_L) for all the 
excess particles to exactly move in phase from one chain to the next 0. In contrast, when 
interaction between particles is switched on we expect that P{t) never reaches 1 for finite 
non-zero r although an oscillatory behavior can still occur in the case of coherent transverse 
hopping. 

Since this issue cannot be adressed by perturbative methods exact diagonalizations of 
small ladder clusters have been performed. Using Lanczos technique, it is easier to calculate 
first the Fourier transform of A{T) defined by 



which can be obtained through a continued fraction expansion. This function is itself quite 
useful because coherence shows up as symmetric peaks around E = 0, while an incoherent 
system is expected to have a broad peak centered around 0. We have calculated these 
quantities for a system of 2 chains with L = 16 sites each, and for a total number of 
particles A^e = 8. For clarity, most of the results presented below correspond to AN^ = 8, 
i.e. to an initial state having all the particles on the same chain, but the main conclusions 
of this work do not depend dramatically on the value of ANe as discussed later. Note that 
antiperiodic boundary conditions have been used in the chains direction to ensure closed 
shell fillings of the corresponding non-interacting systems. Numerical data obtained for 
t±/t = 0.2, ANg = 8 and for various values of J/t ranging from to 2.4 are presented in 
Figures 1, 2 and 3. 

The first important result is that the behaviour depends dramatically on the value of 
J/t. Two typical cases have emerged: 

i) A{E) is a smooth distribution centered around E = 0, and P(r) decreases monoton- 
ically with t. This is for instance what we have obtained for J/t = 0.25 (see Figs, la and 
2b). In this case, the probability that the system goes back to its initial state is continuously 
decreasing, which corresponds to a totally incoherent dynamics. 

ii) A{E) has three narrow peaks, one a.t E = 0, the other two symmetric around E = 0, 
and -P(t) exhibits damped oscillations. This situation is best exemplified by J/t = 2 (see 
Figs, lb and 3). In that case, the probability for the system to go back to its initial state 
reaches again substantial values after going down to almost 0, at least for moderate values 
of T, and the short time dynamics is coherent. Very regular oscillations have also been 
observed for J/t = 0, with a smaller amplitude than for J/t = 2 though (see Fig. 2b). 

The basic properties of a single chain described by the t — J model are known from 
the work of Ogata et al [§. Since we have 8 particles on = 2 x 16 sites, the relevant 
band-filling is n=l/8 (4 particles on each chain). Then, going for J/t = to J/t = 2, the 
exponent a decreases monotonically from 1/8 to [^. Beyond that point, a gap opens in 
the spin sector, although the exact location of the critical value is not known. Finally, a 
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phase separation occurs at Jjt c:^ 2.7. Comparing this behavior with our results (coherence 
for J/t = and 2, incoherence for J/t = 0.25), we conclude that the coherence we have 
detected cannot be related in any simple way to the exponent a. 

We finish the discussion of the data by a comment on the actual role of the parameter 
AA^e- So far we have considered the case ANg/Nf, = 1 where the initial state at r = 
corresponds somehow to a macroscopic perturbation of the system (with full Hamiltonian 
H) far from its equilibrium state. Fig. 4 shows a comparative study of P{t) at J/t = 0.25 
for ANe = 8 and AA^^e = 2. Although ANe = 2 corresponds to a much smaller deviation 
from the absolute ground state (with an excitation energy of only ~ 0.7t compared to ~ 4.3t 
for AA"e = 8) it nevertheless gives rise to a qualitatively similar behavior for P(t) showing 
again no coherence in the transverse motion. 

A natural question which arises then is what causes the special behavior observed only 
for J/t = and J/t = 2. The answer we propose is that these are the only two points 
for which the t — J model is integrable. For J/t = 0, the model is equivalent to the finite 
energy sector of the infinite U Hubbard model, which is known to be soluble by Bethe ansatz 
since the work of Lieb and Wu |T^, while for J/t = 2 an additional supersymmetry again 



makes the model integrable by Bethe ansatz |T1||. The remarkable differences in Fig. 3 
between the results for J/t = 2 and those obtained for J/t = 1.6 and J/t = 2.4 support this 
point of view. For both non integrable cases, the oscillations are much smaller in amplitude 
and involve several frequencies, a sign of a much less coherent dynamics. Coherence might 
be completely lost for arbitrary deviation from the supersymmetric case provided that the 
system size is large enough. 

So the main conclusion of this work is that the short-time dynamics of hopping of elec- 
trons between chains is remarkably coherent if the model describing the chains is integrable 
and more or less incoherent otherwise. 

Let us now compare our results with those obtained previously by other authors. The 
main conclusion of Clarke et al [Q, namely that incoherence can be achieved for a < 1, 
is confirmed by our results for J/t = 0.25, which corresponds to a ^ 0.1. Note that the 



necessary condition t± ^ 2tt—^{vc — Vg) for the observation of coherence proposed by 
Clarke et al is clearly satisfied for J/t = 0.25 because the charge and spin velocities Vc 
and Vg are nearly equal for that particular value of J/t. However, as already been pointed 
out, our results show that coherence does not seem to be linked with the value of a. More 
importantly it seems also that the analysis of the term in the expansion of P{t) for short 
times is not sufficient to study coherence. In Fig. 2a, we have depicted the short time 
behavior of -P(t) for 3 extreme cases: U = (totally coherent), J/t = (partially coherent) 
and J/t = 0.25 (totally incoherent). The curvatures at r = are indistinguishable. 

The main conclusion of Chakravarty and Rudnick namely that incoherence is the 
generic behavior, and coherence is somehow accidental, is also in agreement with our results. 
To go beyond this general statement, one should understand what integrability means in 
terms of the models they have studied, and this is not clear yet. 

Finally the reason why integrability and coherence are related remains an open issue. 



Integrability is already known to have dramatic consequences on the level statistics , the 
distribution being Poisson instead of Wigner in that case, and on the finite temperature 
conductivity |T3| . While the first effect can be qualitatively understood in terms of level re- 
pulsion, a good explanation of the second one is also missing. More work is clearly needed to 
understand the inffuence of integrability on the dynamical properties of correlated electrons. 
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FIGURES 

FIG. 1. Spectral function A{E) for a 2 x 16 t-J ladder with Ne = ANe = 8. Energies and 
parameters are measured in unit of t (i.e. t' = t±/t). (a) Non-integrable case J = 0.25 (and 
J_L = 0.01); (b) Supersymmetric case J = 2. 

FIG. 2. Probability P{t) vs time r calculated at J/t = and J/t = 0.25 for ANe = 8 and a 
transverse coupling t' = t±/t = 0.2. The non-interacting case is also shown (thin dashed line) for 
comparison. Time is measured in unit of the inverse hopping integral 1/t. (a) Small time region; 
(b) Enlarged time interval. 

FIG. 3. Probability P{t) vs time r calculated at the supersymmetric point J/t = 2 and in its 
vicinity for ANe = 8. The non-interacting case is also shown and time is measured in unit of 1/t. 

FIG. 4. Probability P(r) vs time r calculated at J/t = 0.25 and for both ANe = 2 (full line) 
and ANe = 8 (dashed line). 
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